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I. INTRODUCTION 



It was first noticed by Linde in Ref. that broken symmetries may be restored at high 
temperature. Indeed, at very high temperatures, the leading-order corrections are enough to 
restore the broken symmetries in most of the conventional models 0,0]. However, Weinberg 
has shown that, in some models such as the 0(Ni)XO(N 2 ) scalar theory, the symmetry 
is not restored as the temperature increases when the negative coupling between the fields 
extrapolates some limiting value ||. This intriguing phenomenon is called symmetry non- 
restoration (SNR) [or inverse symmetry breaking (ISB)] and has not only academic interest. 
Symmetry nonrestoration at high temperature has crucial consequences in the cosmological 
scenario and could, in principle, avoid the monopole problem Some examples of SNR 

in nature and experimental discoveries may be found, e.g., in |J. Nevertheless, even in 
theories where SNR is believed to be found, there has been some controversy if the phe- 
nomenon, in fact, happens as it was formulated. See, for instance, Ref. |7j for a discussion. 
Recent works show that higher order terms or a nonperturbative study tend to diminish the 
parameter space where symmetry nonrestoration can happen Reference || is a one 

loop self consistent resummation calculation whereas Ref. |7[ takes full account of two loop 
diagrams in a nonperturbative fashion. The parameter space for ISR predicted by the latter 
is bigger than the one predicted by the former. Also the two loop calculation reveals a new 
possibility for ISB where the last term in the breaking sequence is 0(Ni — l)XO(N2 — 1) 
which cannot be seen at the one loop level @. 

In this paper we will apply a nonperturbative method, modified self-consistent resum- 
mation (MSCR), which we have developed recently in order to reinvestigate some 
aspects of SNR. To motivate, let us note that one way to justify the study of resummation 
programs in systems which exhibit symmetry (non)restoration is that those methods repre- 
sent a practical way to recover the reliability of perturbation theory at finite temperatures. 
However, another important reason is that these methods can be a natural manner to take 
into account the medium-modified masses in the computation of several physical processes 
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such as dispersion relations, damping rates and decay widths. The goal of MSCR is to make 
renormalization possible when obtaining the dressed masses in a nonperturbative fashion. 
We will follow a different direction from the conventional approaches about symmetry non- 
restoration at high temperature. These are, in general, concerned mainly with the behavior 
of the modulus of the coupling between the fields in the region of symmetry nonrestoration. 
Various nonperturbative methods have been used to study symmetry nonrestoration, but 
the important issue of renormalization in these kinds of multifield theories has not been dis- 
cussed, as in Ref. ||. Although very important, the details concerning the implementation 
of an appropriate renormalization program within those methods have been overlooked in 
most applications. We focus on the necessary conditions for a resummed multifield theory 
that has its symmetry restored or not at high temperature, to be satisfactorily renormalized. 

To assure renormalizability in a finite temperature application of resummation programs 
to theories which have fields with different masses, one must require that the masses obey 
two essential features. 



(1) A combination that relates them (which can be imposed by some symmetry as in 
the linear sigma model where M% = M% + 2Xu 2 ) JJ. However, as we shall see soon, this 
does not necessarily mean that the symmetry is restored at high temperature. 



(2) The masses in the internal lines of the diagrams have to be the same as in the coun- 
terterms. These requirements are because the resummation solves the infrared divergences 
problem (or the problem of higher order terms in higher order calculations) but may bring 
another one: the divergent temperature dependent parts of the diagrams do not match the 
same terms in the original Lagrangian. That is, the usual constant counterterms cannot 



eliminate terms which are functions of the temperature. At zero temperature [12] or in the 
naive one-loop calculations at finite temperature [[OJ this point does not represent a problem 
since the masses running in the loops are also constants. Then, in resummation methods as 
MSCR it is crucial that the Lagrangian, including the counterterms, suffer the effects of the 



resummation for the model to be satisfactorily renormalizable [[|,|ll|,[14]] . The second neces- 
sary condition is naturally fulfilled by MSCR since the method is based on the recalculation 
of the self-energy using in each step the masses obtained in the previous one. In this way, 
the absorption of the unwanted ultraviolet divergences is always guaranteed once the masses 
running in the loops are necessarily the same as in the counterterms. Regarding the first 
condition, in Sec. Ill we show that we can construct a renormalizable theory (in the sense 
we are discussing) that exhibits both symmetry restoration and nonrestoration, depending 
on the values given to the parameters which govern the relation between the masses. 

This paper is organized as follows. The application of the MSCR method in the 
0(Ni)XO(N 2 ) model is addressed in Sec. II. In Sec. Ill we calculate the dressed masses 
and determine the necessary conditions for a resummed multifield theory which displays 
symmetry restoration or nonrestoration at high temperature to be renormalized. In Sec. IV 
we discuss some numerical results. Section V is devoted to the conclusions. 



II. THE 0(N 1 )XO{N 2 ) MODEL 



Consider the 0(Ni)XO(N 2 ) model described by the Lagrangian 



£ = £ u (0) + £ mt (0) + £ ct (</>) 



(1) 



where 



£ = ^E[(^) 2 --^ 
z i=i 



the interaction Lagrangian is expressed as 



int 



E 

i=l 



A 



4! 



A 

T 



(2) 



(3) 



and the counterterm Lagrangian which is necessary to render the theory finite up to a given 
order is expressed as 



i=i 



(4) 
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The conditions required for the model described by Eq. (|T|) to be bounded from below 

are 

Ai > 0, A 2 > and X 1 X 2 > 9A 2 . (5) 
Let us now apply some topics of MSCR in the 0(Ni)XO(N2) model. The MSCR recipe 



TTfl dictates that 

Ml n = (B hn + l)M^ n _ 1 + Tl t (K , \K\ = 0, M^O, (6) 

where n is the order of the nonperturbative correction, B n is the coefficient of the appropriate 
counterterm and TI^Kq, \K\ = 0, M$ n _i) is the self-energy of the field i in a given order in 
the perturbative expansion at zero external three momentum. This assures the cancellation 
of the ultraviolet divergences since the masses running in loops are necessarily the same 
as in the counterterms. The worse divergences i.e., the infrared, will be cutoff by the 
recalculation (resummation) of the self-energy, which is achieved for n > 1. Then, the 
nonperturbative correction to the thermal mass of the field fa (i — 1,2) to one- loop order, 
which is independent of the external momentum, is given by 



Ml n {T) = (B lin + l)Ml n _ x - {-^-) Kn-i(T)~ (7) 

A N 2 Ml^^l + Finite + ^(Ml^T), Ml n _ x {T)), 



2(4tt 



Ml(T) = (B 2 , n + l)Ml n _ x - ^ (^±^) Ml^T)^ (8) 
A : N 1 M*(T)l + Finite, + nf (M^^T), M*(T)), 



2(4vr) 

where Finitei is the finite part of the divergent contributions, Ilf is the temperature- 
dependent contribution and i = — 7 + log(47r) (this is the modified minimal subtraction 
scheme used in dimensional regularization that we have adopted here), where 7 is the Euler 
constant. 

Let us stop the calculations in the 0(Ni)XO(N2) theory for a while and turn the at- 
tention to a more simple case, in order to compare some results of the MSCR method with 



another approach [plj. In Ref. [|T5 1 , Altherr computed the self-energy up to two-loops in 



the A</> 4 model, and after had identified the origin of the leading infrared part (which came 
from the double scope diagram), he proposed a possible cure for this infrared singularity. He 
suggested as a cure the thermal mass obtained in the massless limit of the one-loop order, 
which could be used as a cutoff. This cutoff mass was put only in the infrared part of the 
self-energy cited above. He saw that this (part of the two-loop order in the perturbative 
expansion) manipulated contribution could be gotten by N-loop diagrams. So he gave up 
the perturbative expansion and considered the most infrared singular diagrams, which are 
the daisy types, as a nonperturbative contribution to the mass shift. He found that the 
result of the summation of these N-loop diagrams is not very different from the perturbative 
approach. 

In Ref. we also treat the A0 4 model up to two loops, but using MSCR. As explained 



before, this method is based on consistently obtaining the pole of the thermal corrected 
propagator, up to a given number of loops in the perturbative expansion, totally free from 
UV and IR divergences, due to the resummation nature of the method. This is completely 
different from what Altherr has done, since with MSCR one does not have to replace the 
perturbative expansion by a nonperturbative approach. Besides, a nonperturbative approx- 
imation has to be done exactly to recover the validity of the perturbation series. As can 
be seen in Ref. [|ll|] , the corrections to the free propagator are the one-and two-loop graphs 
from the perturbative expansion plus infinitely many others obtained nonperturbatively, like 
superdaisy types. MSCR has the advantage of keeping the same fundamental theory when 
calculating the pole of the corrected propagator, using the same counterterms as at zero 
temperature. The only difference is that the mass running in the loops changes in each re- 
calculation (each summation of an infinity set of diagrams) of the self-energy. MSCR is thus 
more natural and keeps the loop expansion, as in Ref. [Q. With the MSCR, the "starting" 
effective Lagrangian can be chosen such that the mass running in the loops is the one-loop 
thermal mass (or, in other words, the plasmon mass m 2 = ^AT 2 ), avoiding in this way 
infrared problems and respecting the symmetries of the Lagrangian, since a loop expansion 



is an expansion in powers of the Lagrangian. 

Then, MSCR can be understood as an alternative and consistent resummation method, 
whose renormalization is in the same way as in the zero temperature case. More simple 
and less complicated naive approaches have found several difficulties when trying to find 
appropriate renormalization conditions for the divergent gap equations, as discussed in ||. 

III. THE DRESSED MASSES IN THE 0(N 1 )XO(N 2 ) MODEL 

We shall now compute the dressed masses and investigate the basic conditions in order 
for the multifield ( "multimasses" ) model to be renormalizable, i.e., completely free from 
"temperature dependent" ultraviolet divergences. To be more clear, (B,i, n + l)M? n _ 1 in 
Eqs. (0) and (|) should be written as {B i>n + l)m? + (B i>n + l)nf en (M 2 n _ 2 ), since in these 
equations Mf n _ x = m 2 + nf en (M 2 n _ 2 ). See the appendix for details. The coefficients of 
the temperature dependent mass counterterms Bi^ n are fixed in a manner to cancel not 
only divergences proportional to ITf en (M 2 n _ 2 ) but also these terms together. This is to 
avoid overcounting of diagrams || . In order for these temperature-dependent infinities to be 
canceled properly by C ct let us consider the most general case of (linear) combination between 
the squared masses. Namely, Mf n „ 1 (T) = aMf^-^T) +J 7 , where a may be a constant and 
T is some well behaved function which can, in principle, depend on (Aj, A, N\,N2,T) and 
may be related with the vacuum expectation value (VEV) of a scalar field. For simplicity 
we will take T as being zero and, as we are going to see below, a particular choice for a 
is a = ±1 with a = 1 representing symmetry restoration and a = —1, on the contrary, 



1 Remembering again the linear sigma model, M 2 = M 2 + J-, with T = 2Az^ 2 (y is the VEV of 
the sigma field), where this relation holds at tree level and in the low temperature region. At 
high temperature J- = 0, signalizing that the symmetry is restored. In the region of intermediate 
temperatures, i.e., around T c (the critical temperature) the one-loop analysis is not enough to 
describe the behavior of the system and the validity of this relation ||. 
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representing symmetry non- restoration. In fact, we will prove in the next section that, at 
sufficiently high temperature and for realistic values of Ni and N 2 and appropriate values 
the couplings respecting the boundness condition, these exact values for a can be found. We 
argue that in ordinary field theory this relation is always satisfied, since at high temperature 
the fields have a T 2 mass behavior. To the best of our knowledge, the exception is left to a 
supersymmetric theory with flat directions which do not have interactions strong enough to 
be in thermal equilibrium at high temperature. A direct consequence of this is the absence 
of T 2 mass terms ||16|| . 

Since n —>■ oo (which makes sense within a resummation program) then n ~ n — 1. 
This implies that (T) = ±M 2 n (T) = ±M 2 (T). After this, the coefficients of the mass 
counterterms at one-loop order are easily found to be 



Let us now see the consequences of this "prescription" on the masses resummed by an 
infinity set of "daisy" and "superdaisy" diagrams which are achieved in the n — > oo limit 



A. First case: Symmetry restoration (a = 1) 

This is the common case where the symmetry is restored at high temperature and will 
simply furnish a relation between the couplings 





fl9yTTH. The study will be divided into two cases. 




(10) 



which at high temperature gives for both positive and negative A, 



S 



Mi + 2) V Nx + 2 

Equation (|ll|) shows the result that for N\ = N 2 we have Ai = A2. Although A here 
can assume negative values, it is not completely free to be larger than the coefficient of T 2 



in Eq. (|T0|). This would imply a nonphysical theory since the third condition imposed for 
the model to be bounded from below in Eq. (El) would not be satisfied, i.e., one would have 



AiA 2 < 9A 2 



N3N2 



(7V!+2)(iV 2 +2)- 



B. Second case: Symmetry nonrestoration (a = —1) 

In this situation we have Mf = — Mf . One can choose Mf as being the mass which has 
the negative coefficient such that 

A 2 ^^ + AiVx < 0, (12) 
so, inverse symmetry breaking will happen for 

W > # m ■ (13) 



Nt V 3 

As we will see below, the exact constraints a = 1 as well as a = —1 can be found for 
realistic values of the parameters of the model, respecting the boundness condition expressed 
by the inequalities @. 

IV. NUMERICAL ANALYSIS 

With the intention of verifying that the approximations imposed on the renormalized 
theory indeed describe very well the cases of symmetry restoration and nonrestoration at 
high temperature let us define new variables as 



Ml = Ml - ml 



2 



r 

24' < i4 > 



Mj = Ml - m\ 



A, 



iV 2 + 2 



+ \N X 



24 



(15) 



The first case, M\ = Mf , that is related to usual symmetry restoration at high tempera- 
ture takes place if one chooses, for example, Ni = 90, iV 2 = 24 (where we have followed Ref. 
H by choosing these values in order the model mimics the Higgs sector of a £77(5) grand 
unified model), A 2 = 0.83 and A = +0.09 which after put in Eq. (O) gives Ai = 0.43 and 
the exact value a = 1, as expected. 

To achieve the second case, we have chosen A2 
We also chose iVi = 90 and N 2 = 24 and Ai = 0.1, A 2 = 0.83 and A = -0.09 such that the 
conditions (H) and (|i~3"D are still satisfied even with the cross coupling A negative. With these 
values chosen, we obtain exactly the wanted result that describes symmetry nonrestoration 



n 2 +2 
3 



I A|iVi < 0, so that Mf is negative. 



iVi+2 



\X\N 2 ] 



a = f — jv 9 +a ; ; — t = — 1 anc ^ consequently, M? = — Mf . In Fig. 1 one sees the behavior of the 

|A 2 — 5 |A|iViJ 

symmetrically opposite masses M\ (> 0) and M| (< 0) showing symmetry nonrestoration 
at high temperature. 

It is important to point out here that: 
(1) typical values for ml and m\ do not alter the results and conclusions since the high 
temperature limit is defined to be T » m\ 2 ; 



M 2 

(2) the model is renormalizable for any value of a = 



with the 



couplings satisfying the inequalities to preserve the boundness condition. In this case, 
the coefficients of the mass counterterms would read 



B2,n = 

as shown in the Appendix. 



Ax ( Ni + 2 
3 

iV 2 + 2 



2(4tt) 2 

A 2 
2(4vr) 2 



+ 



A 



2(4tt) 2 
A 

2(4tt) 1 



No 



1 



Ni a 



1 

e 
1 

e 



(16) 



V. CONCLUSIONS 



In this paper we have obtained the dressed masses of the neutral multifield 0(^1)10(^2) 
model at finite temperature up to one-loop order. We have used MSCR to take into 
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account higher loops contributions and to study in what circumstances the symmetry 
(non)restoration can occur in a renormalized resummed theory. We have found that in 
the case of a multifield theory there is the necessity of the masses to obey a combination 
that relates them. Fortunately the T 2 mass behavior is a natural result of ordinary field 
theory at high temperature which guarantees this constraint. 

We have shown that renormalizable resummed multifield theories which exhibit the par- 
ticular cases of equal masses (symmetry restoration) or symmetrically opposite masses be- 
havior (symmetry nonrestoration) at high temperature can be constructed, depending on 
the physical values given to couplings which govern the strength of the interactions. In the 
future we want to use MSCR beyond the leading order in the perturbative expansion in 
order to verify the existence of other channels of ISB [like 0(Ni — l)XO(N 2 — 1)], as found 
in 0. 
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APPENDIX A: DETERMINATION OF THE COEFFICIENTS OF THE MASS 

COUNTERTERMS 

With MSCR the coefficients of the temperature dependent mass counterterms are deter- 
mined as follows 

M 2 n = (B i>n + l)M? n _ 1 + M^ln-i) (Al) 
(B hn + IK 2 + (B i>n + l)nf e "(M 2 n _ 2 ) + ^(M^O = 

M l,n-1 

(B i>n + l)m 2 + (B hn + l)Ilf n (M 2 „_ 2 ) + a{m\ + nf e ™(M 2 n _J]^ 
+d i \ml + uf n (Ml n _ 2 )]l + Fintte, + flf (M^), 
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where c x = (^3^) an d d i = ~2(tr N ^ for M ln and c 2 = ~2{tf N ^ and d ^ = 

N 2 „ +2 ) for M| n . As can be seen in the equation above, all divergences can be 



2(4tt) 2 V 3 

absorbed, provided M^ n _ 1 (T) = aM% ^(T) + JF. As explained earlier, for simplicity T is 
taken as zero. This implies 

m\ + n«™(Ml 2;n _ 2 ) = \m\ + inf en (M 1 2 2;n _ 2 ), (A2) 

which allows one to get the coefficients of the mass counterterms which read 

B x , n = - (ci + di (A3) 

#L,n = #l,n - 1, V Tl > 1, 

and 

5 2 ,n = - (a c 2 + rf 2 ) i, (A4) 

#2,n = ^2,n - 1, V n > 1. 

It is worth noting here that at zero temperature or in the naive finite temperature calcu- 
lations, there is no necessity of this constraint between the resummed thermal masses since 
the masses running in the loops are constants. This constraint is an effect of the resumma- 
tion approach which consistently forces the masses running in the loops to be temperature 
dependent to fight against the breakdown of perturbative expansion. The renormalization 
of the resummed thermal masses at the end results in that no new counterterms are required 
and the Feynman rules are unchanged. 
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Figure Caption 

FIG. 1: The behavior of Mf. The solid line is Mf and the dashed line is Mf = -M\. 
This result was obtained for realistic values of the parameters of the model. 
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